Engineering the magnetic properties of the Mni3 cluster by doping 
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With a goal to produce giant magnetic moment in Mni3 cluster which will be useful for practical 
applications, we have considered the structure and magnetic properties of pure Mni3 cluster and 
substitutionally doped it with X = Ti, V, Cr, Fe, Co, Ni atom to produce Mni2X clusters. We find 
that Ti and V substitutions in Mni3 cluster are the most promising in terms of gaining substantial 
binding energy as well as achieving higher magnetic moment through ferromagnetic alignment of 
atom-centered magnetic moments. This has been demonstrated in terms of energetics and electronic 
properties of the clusters. For comparison, we have also studied the effect of N-capping of Mni3 
cluster, predicted in the earlier work [Phys. Rev. Lett. 89, 185504 (2002)] as a means to produce 
stable giant magnetic moment in Mn clusters upto cluster size of 5 Mn atoms. 

PACS numbers: 36.40.Cg, 73.22.-f,71.15.Mb 



I. INTRODUCTION 



Behavior of Manganese is very special both in the bulk 
and cluster forms. The electronic configuration of the 
atom has a half filled 3d shell and filled 4s shell. Struc- 
tural and magnetic properties of manganese clusters have 
been studied in the past»i~— It has been shownir— that 
small Mii-clusters containing up to 4-5 atoms exhibit fer- 
romagnetic coupling of individual atomic magnetic mo- 
ments with a magnetic moment per atom comparable 
with the magnetic moment of Mn at the atomic limit, 
which is 5 fiB- Interestingly, the coupling becomes ferri- 
magnetic as the cluster size increases and the net mag- 
netic moment falls abruptly as shown in the experimen- 
tal works^ and also supported by the first principles 
studies*^ The atom centered individual magnetic mo- 
ment, though, remains of the order of atomic moment. 
Therefore, it will be of great technological interest if the 
fcrrimagnetic coupling can be transformed into ferromag- 
netic, thereby resulting into a giant magnetic moment. 
One way to accomplish this is by doping.— Furthermore, 
if the doping also increases the stability of the cluster 
compared to the pure Mn-cluster of same size, then it 
will be an additional achievement on top of spin engi- 
neering. Such clusters will be promising building block 
for making new cluster assembled materials. Along the 
same line of thought, Rao and Jena^ previously studied 
the structure and magnetism of N-cappcd Mn clusters. 
Based on their study carried out on Mn clusters of size 
upto 5 Mn atoms, they proposed that N-capping can be 
used as a viable means to enhance both the stability and 
the ferromagnetic coupling among Mn atoms in a Mn- 
cluster. In the present work, we have proposed an al- 
ternative route and studied systematically the effect of 
substitutional doping of Mni 3 cluster to produce Mn i2 X 
clusters with X = Ti, V, Cr, Fc, Co, Ni atom. We find 



that Ti and V substitutions are the most effective as both 
the magnetic moment and stability increase significantly 
for the Mni2Ti and M1112V clusters compared to the pure 
Mni3 cluster. To compare our finding with the proposal 
of Rao and Jena^ through N-capping, we have also stud- 
ied the structure and electronic properties of Mni 3 N clus- 
ter. Although, the N-capping enhances the stability and 
magnetic moment compared to the pure Mni 3 cluster, 
the gain in total magnetic moment is found to be smaller 
compared to that of Mn^Ti and Mn^V clusters. Also 
the substitutional Ti-doping and V-doping are more in- 
teresting since they keep the volume of the cluster intact 
as opposed to N-capping which increases the volume of 
the original clusters. The rest of the paper is organized 
in the following manner. Section II describes the com- 
putational details. Section III summarizes the results 
obtained in our study. This section consists of five sub- 
sections. The sub-section A is devoted to the results de- 
scribing pure M1113 cluster. The sub-sections B and C are 
devoted to results on M1112X clusters. In sub-section D, 
we analyze the origin of enhanced stability and magnetic 
moments in Mn^X clusters. Sub-section E describes our 
comparative study on N capped Mn 13 cluster. We end 
the paper with a conclusion in section IV. 



II. 



COMPUTATIONAL DETAILS 



Our calculations have been performed using density 
functional theory, within the pseudo-potential plane wave 
method as implemented in the Vienna ab initio sim- 
ulation packaged We have used the projector aug- 
mented wave pseudo-potential a 14 ' 1 5 and Perdew-Bruke- 
Ernzerhof exchange-correlation functional^ for spin- 
polarized generalized gradient approximation. The wave 
functions were expanded in the plane wave basis set with 



the kinetic energy cut-off of 270 eV and reciprocal space 
integrations were carried out at the T point. Symmetry 
unrestricted geometry optimization has been performed 
using the conjugate gradient and the quasi-Newtonian 
methods until all the force components became less than 
a threshold value of 0.005 eV/A. Though the noncollinear 
effect of atomic spins is significant for Mn-clusters, it has 
been shown in the earlier workii that for Mn 13 cluster, 
the degree of noncollincarity is very small. Therefore, we 
have optimized the structures for all the spin multiplic- 
ities M = 2S + 1 under the approximation of collinear 
atomic spins arrangement to determine the magnetic mo- 
ment of the minimum energy structure. Simple cubic 
super-cells were used with the periodic boundary condi- 
tions, where two neighboring clusters were kept separated 
by at least 12 A vacuum space to make the interaction 
between the cluster images negligible. The binding en- 
ergy (Eb) of each Mn^X cluster has been calculated with 
respect to the free atoms as 

E B (Mn 12 X) = 12£(Mn) + E(X) - E(Mn 12 X) (1) 

where £?(Mni 2 X), E(Mn) and E(X) are the total energies 
of the Mni2X cluster, an isolated Mn atom and an iso- 
lated X atom respectively. We also have analyzed the gap 
between the highest occupied molecular orbital (HOMO) 
and the lowest unoccupied molecular orbital (LUMO) for 
the magnetic clusters by defining the spin gap as follows: 
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The system is considered to be stable if both Si and 5 2 
are positive. Detailed discussion about 5i and S 2 can be 
found elsewhere *2r— 



III. RESULTS AND DISCUSSIONS 
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FIG. 1: (Color online) Plot of total binding energies for icosa- 
hedral, hexagonal closed pack and cuboctahedral structures 
of the Mni3 cluster as a function of total magnetic moment 
(M) which is the difference between the number of electrons in 
the up (i.e Nf) and down (i.e NJ,) spin components. The data 
point corresponding to the minimum energy Mni3 cluster has 
been encircled. 




FIG. 2: (Color online) Atomic spin configurations for the min- 
imum energy structure of pure Mni3 cluster. The cyan (light 
gray) colored balls correspond to atoms with positive mag- 
netic moment and blue (dark gray) colored balls correspond 
to atoms with negative magnetic moment. The convention 
for numbering the atoms is shown, which is used in Fig. 7. 



A. Pure Mni3 cluster 

Since the transition metal clusters generally prefer 
compact geometry to maximize the interaction between 
the rather localized d — d orbitals,— we have considered 
icosahedral, cube-octahedral and hexagonal closed pack 
(HCP) structures as the most probable starting geome- 
tries to determine the minimum energy structure (MES) 
for Mni3 cluster. These three are the most compact 
and highly coordinated structures for the 13 atoms clus- 
ter. Each structure was allowed to relax for all possible 
collinear spin configurations of the atoms. Fig.l shows 
the variation of total binding energy as a function of to- 
tal magnetic moments for these three structures of Mni3 
cluster. 



We found an icosahedral structure with total magnetic 
moment of 3 fiB and total binding energy of 28.023 eV is 
the MES for the Mni3 cluster, as predicted earlier J£ An 
icosahedron of 13 atoms can be viewed as consisting of 
two pentagonal rings connected with each other through 
the central atom and capped with two apex atoms at its 
ends as shown in Fig. 2. We also observed a spin seg- 
regated state of atomic moments in the MES of Mni3 
cluster in the sense that the atoms within each of the 
two pentagonal rings of icosahedral Mni3 cluster are fer- 
romagnetically aligned, while the two pentagonal rings 
are anti-ferromagnetically coupled with each other (cf. 
Fig. 2). This is in agreement with the previous report by 
another group J- 
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FIG. 3: (Color online) Plot of total binding energies of icosahedral Mni2X clusters as a function of total magnetic moment (M). 
The total binding energies are plotted with respect to the binding energy of the MES of icosahedral Mni3 cluster. Panels (a) 
and (b) correspond to the substitutions by elements left and right of Mn in the periodic table respectively. The dotted lines 
through zero in panels (a) and (b) correspond to the binding energy of MES for pure M1113 cluster. The data points around 
the minimum energy structure of each Mni2X cluster have been emphasized in the insets with the MES by specially marked. 



B. Energetics of Mni2X clusters 

For the calculations reported in the following, we have 
considered doped Mni 2 X clusters of the icosahedral sym- 
metry with the impurity X atom sitting at the center of 
icosahedron and optimized it for all possible collinear spin 
configurations of atoms. We have considered six differ- 
ent substitutions at the central Mn site of Mni3 cluster: 
considering the substituting X atom as one of the three 
elements left of Mn in the periodic- table (Ti. V, Cr) and 
of the three elements right of Mn in the periodic table 
(Fe, Co, Ni). This assumed geometry of icosahedral sym- 
metry with the impurity X atom sitting at the center of 
icosahedron needs some discussion. In our earlier work on 
V-doped Co clusters^ we have also shown by ab-initio 
molecular dynamic study that the ground state structure 
of C012V cluster, adopts an icosahedral pattern with the 
substituting V-atom always preferring to sit at the cen- 
tral site of icosahedral C012V cluster. Theoretical works 
on clusters with isosahedral symmetry with substitution 
at the central site have been reported for MC012 (M=Ti, 
V, Cr, Mn, Fe, Co and Ni) clusters^ Fei 2 X (X = Be, 
O, Mg, Al, Si, S, Ti, Mn, Fe, Ni, Zn, Se, Pd) clusters^ 
and AI12X (X=B, Al, Ga, C, Si, Gc, Ti, As) clusters^. 
Therefore, there exists enough reports in the literature 
considering the substitution at the central site by atoms 
of similar or larger sizes, justifying our choice of the im- 
purity A-atom doped at the central site of icosahedral 
Mni3 cluster. Furthermore, as discussed later in section 
III.D, Mn-X dimcrs were found to have higher binding 
energies compared to Mn2 dimer. This in turn would 
indicate that single atom substitution of Mn by X atom 
will prefer a configuration where the number of Mn-X 
bonds are maximized, which is possible for the substitu- 
tion at the central position. Our total energy calculations 



by putting the substituting atoms at different sites of the 
icosahedral Mni3 cluster support this idea. 

Fig. 3 shows the plot of binding energies as a function 
of total magnetic moments (Nt - NT.) of all the Mn^X 
clusters. In order to see the effect of substitutional dop- 
ing distinctly, we have plotted the binding energy for 
each spin multiplet of Mn^X clusters with respect to 
that of the MES for pure Mni 3 cluster i.e E B (Mni 2 X)- 
Eb(Miii3)m.es- Positive value of this quantity means 
gain in binding energy upon substitutional doping. It 
is seen from Fig. 3(a) that gains in binding energies are 
achieved in case of Ti and V substitutions for a range of 
values of total magnetic moments, from ~ 2 fis to ~50 
fiB, while for Cr-substitution there is no gain in binding 
energy except a small gain for the magnetic moments of 
22 hb and 10 {1b- On the other hand, for Fe, Co and 
Ni substitutions, there is also gain in binding energy for 
each of them upto a total magnetic moment of around 30 
Hb as seen from Fig. 3(b). Notice that these gains are not 
as high as we observed in case of Ti and V substitutions. 
However, they are better than the Cr-substitution. 



C. Energetics and magnetic moments of MESs of 
Mni2X clusters 

Table U shows the total binding energies and total mag- 
netic moments of the MESs of Mn^X clusters. It is 
seen that the MESs of all the doped clusters have higher 
binding energies and magnetic moments as compared to 
those of the MES for Mni3 cluster. We note again that 
the MESs of the Mn i2 Ti and Mm 2 V clusters have sig- 
nificantly higher binding energies than Mni3. Similarly, 
the total magnetic moments are found to be significantly 
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FIG. 4: (a) Binding energies [Ediff = Es(Mni2X)MBS - Eb (M1113 )mss] and (b) magnetic moments [Mdiff = M(Mni2X)A/ES 
- M(Mni3)MBs] of the minimum energy structures of M1112X clusters with respect to those of the minimum energy pure M2113 
cluster respectively. The straight lines through zero in panels (a) and (b) are the reference fixed at the binding energy and 
total magnetic moment respectively of the MES of pure Mni3 cluster. 



TABLE I: Total binging energies and total magnetic moments 
(M) of the minimum energy icosahedral structures of Mni2X 
clusters, X = Ti, V, Cr, Mn, Fe, Co, Ni. 
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large for the MESs of clusters doped with elements to 
the left of Mn in the periodic table. Moreover, there 
are several isomers [as seen from Fig. 3(a)] of Mn^Ti 
and Mni2V clusters with total magnetic moment in the 
range of 3-50 (1b, which have gain in binding energies 
compared to the MES of pure Mni3 cluster. Combin- 
ing the results of Fig. 3 and Table 1, it is therefore seen 
that the substitutional Ti-doping as well as V-doping of 
the Mni3 cluster are the most preferable in view of gain- 
ing significant amount of both the binding energy and 
the total magnetic moment. For the MES of Mn^Cr 
cluster, although the gain in magnetic moment is quite 
significant, the gain in binding energy is the least. Sim- 
ilarly, Mni2Fe, Mn^Co and Mni2Ni clusters also have 
some isomers with total magnetic moments in the range 
of around 3-30 /j,b [as seen from Fig. 3(b)], which have 
both gains in binding energy and magnetic moment with 
respect to those of the MES for pure Mni 3 cluster. How- 
ever, the MESs for Mn^Fe, Mn^Co and Mn i2 Ni clusters 
have total magnetic moments of 8 (1b, 7 [Ib and 8 \xb re- 
spectively, which are higher than that of the MES for 
M1113 cluster too, but significantly lower than that of the 
MESs for Mni2Ti and Mn^V clusters. In order to have 
a visual depiction of the trend in the choice of substitut- 
ing elements in a concise manner, the binding energies 
and magnetic moments of the MESs of Mn^X clusters 



with respect to those of the minimum energy pure Mni3 
cluster respectively i.e E B (Mni 2 X) M _ES - E B (Mni 3 ) M £;s 
and M(Mni 2 X)j\/£;s - M(Mni 3 ) A /,Es are plotted in Fig||| 

In FigEJ we have plotted the total density of states and 
center atom projected density of states of the MES of 
each Mni2X cluster. As is seen from FigJSJ the difference 
between the majority and minority spin states is large 
for the optimal icosahedral Mn^Ti, Mn^V and Mn^Cr 
clusters, which results in large magnetic moments for 
these clusters. On the other hand, the magnetic moments 
for the optimal icosahedral structures of M1113, Mn^Fe, 
M1112C0, Mni2Ni clusters are relatively small, as the ma- 
jority and minority spin states mostly cancel each other. 
Center-atom projected density of states shows that it 
has finite but small magnetic moment in case of MESs 
of Mn i2 Ti, Mni 2 V and Mn i2 Cr clusters. On the other 
hand, its magnetic moment is vanishingly small for the 
MESs of Mni 3 , Mn^Fe, M1112C0 and Mn i2 Ni clusters. 



D. Origin of stability and enhanced magnetism of 
the MESs of Mni 2 X clusters 

To understand the relative stability of the doped 
Mn 12 X clusters, we have computed the binding energies 
of MnX dimers as given in Table II and also plotted 
the spin gaps [as defined in EqJ5] for the MESs of the 
M1112X clusters in Fig[6] Our calculated binding ener- 
gies for dimers follow the same trend as we observed in 
case of clusters. It is seen from Table II that MnTi dimcr 
has the maximum binding energy of 3.14 eV. Again the 
binding energy of MnV dimer (~ 2.72 eV) is also quite 
high. Therefore, the large gains in binding energies of 
Mni2Ti and Mn^V clusters come from the large bind- 
ing energies of MnTi and MnV dimers respectively. On 
the other hand, the MnFe, MnCo and MnNi dimers also 
have large binding energies compared to that of the Mn2 
dimer. Consequently, the Mn^Fe, Mn^Co and Mn^Ni 
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FIG. 5: (Color online) Plot of total density of states (left 
panel) and center atom projected total density of states 
(right panel) of the minimum energy icosahedral structures 
of Mni2X clusters. Black curves correspond to the results 
for majority spin channel and red dotted curves for minority 
spin channel. The dotted line through zero parallel to y-axis 
in each system corresponds to Fermi-energy of that system. 
Gaussian smearing (0.1 eV) has been used. 



clusters also have larger binding energies compared to the 
pure Mni3 cluster. 

For comparison, we have also listed the previously cal- 
culated values of binding energies and bond lengths of 
the dimcrs as reported in Ref.26. Our calculated values 
follow the similar trend as obtained from the reported 
values for both bondlcnghts and binding energies of the 
dimers. However, we note that the reported binding en- 
ergies in Ref. 26 show the least value for MnCr dimer 
and followed by Mn2 dimer, while our calculated binding 
energies for dimers show the least value for Mn2 dimer, 
followed by MnCr dimer. In order to verify our calcu- 
lation for MnCr dimer, we have also performed the cal- 
culation for Cr2 dimer. Our calculated binding energy 
and bond length for Cr2 dimer are 2.08 eV and 1.59 A 
respectively. These values are consistent with the ex- 
perimental values for Cr2 dimeri^i Considering a mixed 
dimer with Mn and Cr, then the corresponding values for 
the MnCr dimer should be in between the corresponding 
values of the two pure dimers - Mn2 and Cr2. Our calcu- 



TABLE II: Binding energies and bond lengths of MnX dimers 
in the present calculation. For comparison, we have also listed 
the corresponding values as reported earlier (Ref. 26) for all 
the dimers. 



Dimers Binding energy (eV) Bond length (A) 
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Ref. 26 
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MnTi 


3.14 


2.76 


1.80 


1.76 


MnV 


2.72 


2.88 


1.65 


1.69 


MnCr 


1.17 


0.93 


2.34 


2.46 


Mn 2 


1.01 


1.15 


2.58 


2.62 


MnFe 


1.76 


1.57 


2.39 


2.42 


MnCo 


2.17 


1.91 


2.06 


2.09 


MnNi 


2.67 


2.41 


2.06 


2.09 



lated values of binding energy and bond length for MnCr 
dimer follow this trend. Both the spin gaps (Si and ^2) 
as plotted in Fig|51 are also positive for all the clusters 
studied here indicating the formation of stable clusters. 
We find that the MESs of Mn 12 Ti and Mn 12 V clusters 
have large positive value of 5's, which in turn, again indi- 
cate the high stability of the Mn^Ti and Mn^V clusters 
compared to the others. 

To understand the magnetic behavior of the MESs 
of the Mni2X clusters, we have calculated the atom- 
centered magnetic moments within a specified sphere 
around each atom of the MESs of doped Mni2X clus- 
ters and pure Mni3 cluster. Fig. [7] shows the pictorial 
representation of our analysis. Focusing first on the in- 
dividual moments, we find that the moments of all the 
atoms for each cluster remain close to the atomic moment 
(~ 3-4 /is) with the exception of the highly coordinated 
central atom (atomic index 13) which has a small magni- 
tude for all the Mn 12 X clusters. It is therefore the align- 
ment of individual moments of the surfa.ee atoms which 
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FIG. 6: Plot of spin gaps [as defined in Eq. 2] of the MESs 
of Mni2X clusters. 
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Atomic indices 

FIG. 7: (Color online) Individual atomic magnetic moments (m) of the MESs of Mni2Ti, Mni2V, Mni2Cr, Mni3, Mni2Fe, 
Mni2Co and Mni2Ni clusters. The atomic indices 1-12 correspond to the surface atoms (same for each Mni2X cluster), while 
index '13' corresponds to the central atom, as shown in Fig. 2. Lengths of the bars correspond to the magnitude of atomic 
moments. Red colored (dark) bars in the positive direction correspond to +ve magnetic moments and green colored (gray) 
bars in the negative direction for -ve magnetic moments. The central site is the most coordinated site and therefore has very 
small magnetic moment for all the clusters. 



plays the major role in determining the total magnetic 
moment of the MESs of the Mn 12 X clusters. In case of 
doping with atom of elements left of Mn in the periodic 
table (i.e for the Mni 2 Ti, Mni 2 V and Mni 2 Cr clusters), 
it is seen that the surface Mn-atoms are mostly ferro- 
magnetically aligned. For example, only the moments of 
two surface atoms are antifcrromagnctically aligned with 
the moments of other surface atoms in case of MESs for 
Mn 12 Ti and Mn 12 V clusters. Similarly, moments of three 
surface Mn atoms are antiferromagnetically aligned with 
the moments of other surface atoms in case of the MES 
for Mn 12 Cr cluster. On the other hand, the surface Mn 
atoms are mostly antiferromagnetically coupled among 
themselves in case of substitution with atom of the ele- 
ments right of Mn in the periodic table (i.e for Mni 2 Fc, 
Mn 12 Co and Mn 12 Ni clusters) and the net magnetic mo- 
ment for each of these clusters is small, close to the total 
magnetic moment of two ferromagnetically coupled Mn 
atoms which is 10 (Xb^ 

We further find that the magnetic moment of the cen- 
tral atom, though it is small compared to that of the sur- 
face atoms for each doped cluster, varies from one doped 
cluster to another. It shows a systematic increase as one 
moves from Mn 12 Ti— >Mni 2 V— s>Mni 2 Cr and attains tiny 
values for Mni3, Mn 12 Fe, Mn 12 Co and Mn 12 Ni. The mo- 
ment at the central site is decided by two factors - (i) 
gradual filling of electronic states in moving from left 
to right of the periodic table and (ii) the hybridization 
with neighboring atoms. The gradual filling of electronic 
states should increase the moment in moving from Ti to 
V to Cr attaining a maximum value for the half-filled 
case of Mn, followed by a gradual decrease in moving 
from Fe to Co to Ni. While the above mentioned trend 
is observed for the central atom of Mn 12 Ti, Mn 12 V and 
Mn 12 Cr clusters, the moments for the central atom are 
all tiny for Mni 3 , Mn i2 Fc, Mni 2 Co and Mni 2 Ni clusters, 
in contrast to the electron filling trend which would de- 
mand the moment to be the maximum for Mni3. This 
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FIG. 8: Plot of centre to surface average bondlengths ((r)) of 
the MESs of Mm 2 X clusters. 



is, however, explained by the second effect, namely the 
hybridization effect. As the bond length decreases, gen- 
erally the magnetic moment also decreases through the 
increase in effective hybridization. In order to examine 
the effect of hybridization, we have, therefore, calculated 
the centre to surface average distance for the MESs of 
all the seven Mn 12 X clusters with X = Ti, V, Cr, Mn, 
Fe, Co, Ni. Fig. [5] shows the plot of our calculation. We 
find the average bond lengths to grow in moving from 
Mn 12 Ti to Mn 12 V to Mni 2 Cr and then dropping sud- 
denly to a smaller value for Mni3 which stays as smaller 
value for Mni 2 Fe, Mn 12 Co and Mni 2 Ni clusters also. The 
hybridization effect therefore takes over the filling effect 
in producing small magnetic moment at the central sites 
for Mn i3 as well as Mn i2 Fe, Mni 2 Co and Mn i2 Ni. There- 
fore, we find that the central atom behaves as magnetic 
atom in case of Mni 2 Ti, Mni 2 V and Mni 2 Cr clusters, 
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while it behaves as essentially nonmagnetic atom in case 
of M1113, Mni 2 Fe, Mni 2 Co and Mni 2 Ni clusters. Con- 
sidering the alignment of moments for both the surface 
atoms as well as central atom, we can then conclude that 
when the central atom is magnetic, the surface Mn-atoms 
are mostly antiferromagnetically coupled with the central 
atom, causing the effective ferromagnetic coupling among 
the surface atoms which gives rise to large total magnetic 
moment of the clusters. On the other hand, when the 
central atom is nonmagnetic, the surface magnetic mo- 
ments interact with each other through super-exchange 
path involving the central nonmagnetic atom, resulting 
into antifcrromagnetic coupling within each other. This 
gives rise to a small net magnetic moment of the whole 
clusters. 



E. N-capped M1113N cluster 

Finally, in order to see the effect of N-capping in M1113 
cluster along the same direction of Ref.12, we have also 
studied the structure, stability and magnetic properties 
of an icosahedral Mn^N cluster with N atom at a trian- 
gular face cap as shown in the inset of FigJ5J To deter- 
mine the magnetic moment of the MES and to identify 
the possible isomers close to MES, we have again opti- 
mized the structure for all possible collinear spin config- 
urations of atoms. We find an icosahedral structure of 
total magnetic moment of 12 /is is the MES for Mn^N 
cluster. We also find that the Mn-N-Mn bond angles of 
N atom with the three nearest Mn-atoms in case of the 
MES of M1113N cluster, are 90.08°, 93.95° and 90.97°. 
This nearly 90° bond angle between Mn-N-Mn has been 
predicted as guiding rule for the ground state structure of 
N-caped Mn-clusters in the earlier work.— It results from 
the interaction of the p-orbitals of N-atom with the s-d 
hybridized orbitals of Mn atoms. As the N atom is added, 
the binding energy of Mn^N cluster improves signifi- 
cantly compared to that of the MES for pure Mni3 clus- 
ter. The energy gained in adding a N atom to M1113 clus- 
ter, is defined as A = -[E(Mn 13 N)~E(Mn 13 )-E(N)]. 

FigJS] shows the variation of A as a function of total 
magnetic moment of the Mn^N cluster. The MES of 
Mni3N cluster has binding energy gain A = 6.42 eV. 
There are also several closely lying isomers with total 
magnetic moments of 22 fis, 10 fJ-B and 2 [Ib which are 
4 meV, 10 meV and 14 meV away from that of the MES 
respectively. So A's for them are also very close to that 
of the MES. These values of A's are consistent with the 
energy-gain predicted in Ref.12 for Mns cluster, which 
is 5.61 cV. In the earlier work,— it has been shown that 
addition of one N-atom to at least five Mn atoms is en- 
ergetically favorable. Here we sec that addition of sin- 
gle N atom is also energetically favorable even for the 
M1113 cluster. This enhancement in bonding due to N- 
addition has been demonstrated earlier— by examining 
the bonding in case of dimcr - Mn 2 versus Mn 2 N. Lack 
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FIG. 9: (Color online) Variation of energy gain A as defined 
in the text as a function of total magnetic moments of the 
Mni3N cluster. The insets show the atomic spin orientation 
of the minimum energy icosahedral structure and the first 
isomer of Mni3N cluster. The cyan (light gray) colored balls 
correspond to atoms with positive magnetic moment and blue 
(dark gray) colored balls correspond to atoms with negative 
magnetic moment The N-atom [represented by the smaller 
sized red (dark gray) colored ball] in case of both the struc- 
tures of inset, have negative magnetic moment. The numbers 
in the parentheses for the insets correspond to total magnetic 
moment and the relative energy to the MES of Mni3N cluster 
respectively. 



of hybridization between the s and d electrons due to 
filled 4s shell of Mn prevents Mn 2 from forming a strong 
bond. As N is attached, the 4s 2 electrons of Mn interact 
with the 2p 3 electrons of N resulting in strong coupling. 
The inset of Fig |H] shows the orientations of atomic spins 
of the minimum energy structure of magnetic moment 
12 /is and the first isomer of magnetic moment 22 fis 
of Mni3N cluster. The local effect of N-capping in now 
very obvious. Effectively, it couples the atomic spins fer- 
romagnetically in the triangular face where it sits as well 
as in the pentagonal ring nearest to it (upper pentagonal 
rings of the two structures in the insets of Fig[9|) . There- 
fore, two pentagons for each of the two structures (MES 
and the closely lying first isomer) in the inset of FigJSl are 
now effectively ferromagnetically coupled. So in addition 
of enhancing binding, N-capping also increases the mag- 
netic moment compared to that of the pure Mni3 cluster. 
However, it is significantly lower than that of the MESs 
for Mni 2 Ti and Mni 2 V clusters, presumably due to the 
fact that N-capping produces a local effect while the Ti 
or V substitutions at the central site effects all the other 
12 surface atoms thereby producing a global effect. 
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IV. SUMMERY AND CONCLUSIONS 



Total magnetic moment of a M1113 cluster is small due 
to anti-ferromagnetic alignment of individual atomic mo- 
ments. Our study shows that it is possible to obtain 
stable ferromagnetic alignment of atomic spins by sub- 
stitutional Ti-doping and V-doping in Mni3 cluster re- 
sulting into giant total magnetic moment. We have com- 
pared our findings with the previously proposed route of 
TV-capping to Mn-clusters. The gain in magnetic mo- 
ment compared to pure Mni3 cluster is found to be more 
in case of Ti and V substitutions. Also substitutional 
doping keeps the volume of the parent cluster intact, 



while the N-capping increases the effective volume of 
the cluster. Therefore, we propose Ti-substitution and 
V-substitution as promising alternative to engineer the 
magnetic structure for Mni3 cluster. 
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